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ABSTRACT 



Recent Hi observations reveal that the discs of spiral galaxies are surrounded 
by extended gaseous haloes. This extra-planar gas reaches large distances (several 
kpc) from the disc and shows peculiar kinematics (low rotation and inflow). We have 
modelled the extra-planar gas as a continuous flow of material from the disc of a 
spiral galaxy into its halo region. The output of our models are pseudo-data cubes 
that can be directly compared to the H I data. We have applied these models to two 
spiral galaxies (NGC 891 and NGC 2403) known to have a substantial amount of extra- 
planar gas. Our models are able to reproduce accurately the vertical distribution of 
extra-planar gas for an energy input corresponding to a small fraction (<4 %) of the 
energy released by supernovae. However they fail in two important aspects: 1) they do 
not reproduce the right gradient in rotation velocity; 2) they predict a general outflow 
of the extra-planar gas, contrary to what is observed. We show that neither of these 
difficulties can be removed if clouds are ionized and invisible at 21 cm as they leave 
the disc but become visible at some point on their orbits. We speculate that these 
failures indicate the need for accreted material from the IGM that could provide the 
low angular momentum and inflow required. 

Key words: galaxies: kinematics and dynamics — galaxies: individual: NGC 891, 
NGC 2403 — galaxies: haloes — galaxies: ISM — galaxies: evolution — ISM: kine- 
matics and dynamics 



1 INTRODUCTION 

The evolution of spiral galaxies is influenced both inter- 
nally by star formation and externally by the environ- 
ment. Whether the former or the latter mechanism is 
dominant is diflicult to determine. Star formation is re- 
sponsible for blowing gas outside the galactic discs, per- 
haps poll uting the Intergalactic Medium (IGM) wit h met- 
als (e.g. iMcCarthv. van Breueel fc HeckmanI Il987ll . Out- 
flows of gas from galactic discs with velocities of the or- 
der of 100kms~^ are observed both in the neutral (e.g. 
iKamphui s. Sancisi fc van der Hulst' 'l99l") and the ionised 
phase (e.g. .Fratern ali, Oosterl oo fc S ancisi 2 004ai) . On the 
other hand, accretion of unpolluted material from the 
IGM is predicted by chemical evolution models to solve 
the so-called G-Dwarf pro blem in the Milky Way (e.g. 
iRocha-Pinto fc Maciei ll996'1. The low metallicities of some 
of the High Velocity Clouds (HVCs) dWakker et al.lll99g^ 
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suggest that cold material is indeed accreted by spiral galax- 
ies like our own llOordll970D . 

In the last years, the study of the connection between 
galactic discs and haloes has flourished due to the collec- 
tion of new high-sensitivity data at different wavelengths. 
Deep Hq observations have revealed the presence of ex- 
tended layers of diffus e ionised gas (DIG) around edge - 
on spiral galaxies (e.g. iHoopes. Walterbos fc Randl 11999) . 
X-ray observations with the new satellites have shown 
the presence of hot gas at distan ces of tens of kpc from 
the plane of galactic discs (e.g. IWang et al.l l200lh . HI 
observations of several spiral galaxies have revealed ex- 
tended th ick layers or haloes surrounding the galacti c discs 
dSwaters. Sancisi fc van dor Huls t 1997.; Matthews fc WoodI 
l2003l) . The kinematics of this extra-planar gas is charac- 
terized by: 1) a d ecrease in rotation velocity i n the 
vertical direction l|Swaters. Sancisi fc van der Hulsd Il997l : 
[Fraternali et al."2005) ; 2) vertical motions from and towards 
the disc^Boonism^^^l,^005J ; 3) a possible general radial 
inflow jFraternah et alJl200ll) 

These studies show that the extra-planar gas is likely 
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to be the result of a complex exchange of gas between the 
disc to the halo regions. In the galactic fountain scheme, 
gas is pushed up by stellar activity, tra vels through the 
halo and eventually falls back t o the disc iShapiro fc Fiek 
ll976l:lBregmanlll980l) . Recentlv. ICollins. Benjamin fc Ra" 
l^2002^ presented a ballistic model of a galactic fountain and 
compared it with Ha long-slit observations of two edge-on 
galaxies: NGC 5775 and NGC891. They found a disagree- 
ment, especially for NGC 891, between the rotational veloc- 
ities above the plane predicted by their model and those 
derived from the data. Other studies ha ve concentrated on 
the modelling of th e multi-phase ISM l)R,osen fc Bregm'aiil 
[1995'; 'Avillcz '2O OOI), the cooling of the hot coronal gas 
CoAell^^^^geteil Il988i ) and the effect of turbulence 
ruck fc Smi 

Some authors have tried to model the extra-planar 
:as as a stationar y medium in hy drostatic equilibrium 
Beniaminll2002ll : iBarnabe et al.ll2005^ . This approach leads 
to solutions that can reproduce the observed gradient 
in the rotation velocity above the plane if the pressure 
of t he mediu m does not depend on the density alone 
iBarn abe ct aO HooH). However, the temperatures of the 
extra-planar (hydrostatic) gas are of the order of 10^ K and 
it is unclear how this medium can be related to cold neutral 
gas. 

Here, we present a new model for the neutral extra- 
planar gas as the result of a continuous galactic fountain 
flow. We use a different strategy from the previous attempts 
mentioned above. We directly compare our dynamical mod- 
els with HI observations (data cubes) by producing, as an 
output of the model, a pseudo-data, cube with the same to- 
tal flux and resolution as the data. The direct comparison 
with the observations is crucial because it removes the in- 
termediate steps of the data analysis such as the derivation 
of rotation curves of the extra-plana r gas or the separation 
of the halo gas from the disc gas (see iFraternali et al.ll20o3) 
that can introduce significant errors. 

We apply this model to the observations of two spiral 
galaxies that currently provide the best examples of detec- 
tion of extra-planar gas: the edge-on galaxy NGC 891 (Sec- 
tion|^ and NGC 2403, which is seen at 63° along the line of 
sight (Section 2J. The study of galaxies viewed at different 
inclination angles is essential to remove projection effects 
and obtain the full three-dimensional picture of the distri- 
bution and kinematics of the extra-planar gas. The data for 
an edge-on galaxy like NGC 891 give information on the ex- 
tent of the gaseous halo and the rotation velocity of the gas 
iSwaters. Sancisi fc van der Hulsdll997i) . In addition, the 
study of a galaxy at intermediate inclination like NGC 2403 
provides information on non -circular (inflow/outflo w) mo- 
tions of the extra-planar gas JPraternah et alJl200ll) . 



2 THE MODEL 

We have modelled the extra-planar gas in a spiral galaxy as 
a continuous flow of material from the disc to the halo. Indi- 
vidual gas clouds are modelled as non-interacting particles 
that are pushed up from the plane of the disc. The trajec- 
tory of each particle is integrated in the galactic potential 
until it falls back to the disc. Runs of the models were also 
performed in which the particles cross the plane once but 



are stopped on their second passage. Positions and veloci- 
ties along these trajectories are projected along the line of 
sight at each time interval, producing a pseudo-cube of sky 
positions and velocities to be compared with the HI data 
cubes. 



2.1 The potential 

The galactic potential and forces were evaluated numerically 
on a grid in the {R,z) meridional plane, with cells 0.1 kpc 
on an edge close to the plane and 0.5 kpc on a side further 
out. The volume densities of both the stellar and the gaseous 
components are given by 



Pdisc(-R, 2) 



Poe 



-R/Rd 



(1) 



where po is the central density, Rd and are the scale 
length and scale height of the disc, and ("(2) is the vertical 
density distribution, either an exponential or a sech^ law. 
We assume that the scale height is independent of R. 

The first part of the r.h.s. of eq. 1 is proportional to the 
surface density of a razor-thin exponential disc. By mod- 
elling the disc as an in finitely fiattened homoeoid, one can 
show JCuddeforJl993l) that its potential is 



<&thm(-R, 2) = 



4G'Eo 1"^. . I 2a \ „ / a 
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S++S- 
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where Eo is the central surface brightness, Kq is the modified 
Bessel function of zeroeth order, and S± = y^z"^ + [ai: R)^. 
The potential of the thick exponential disc (eq. is 



<I>thick(i?,2) = / d2'C(2')$thin(i?,2), (3) 

where <l>thin(J?, 2) is as in eq. |21but with 



The ellipsoidal components (bulge and dark matter 
halo) h ave been mo delled using a double power law density 
profile jPehnen fc Bi nncv 199Si): 

PdAR,z) = poApiy — j \}^~) ' 

where po,dpi is the central density, a is the scale length, 
7 and /3 are the in ner and outer slopes respectively, and 
m = Y^-R^ + z^/q'^ with q the axis ratio, related to the ec- 
centricity e hy q = y/1 — e^. 

The potential and forces of the ellipsoidal compo- 
nents are calculated usin g formulae (2-84b) and (2-88) in 
iBinnev fc Tremaind lll987il with the substitution of the den- 
sity p with that shown in our eq. 01 

The double-power-law profile allows us to model dark 
matter (DM) haloes with different shapes such as pseudo- 
isothermal, 7 = and /3 = 2; NF W profile, 7 = 1 and /3 = 3 
jNavarro. Frenk fc Whit3ll997^ : Hernquist profile, 7 = 1 
and P = 4 llHernauistlll990l) . We have also used these pro- 
files to model bulges that satisfy de Vaucouleurs' i?^''''-law 
(de V aucouleurs 1948). The associated luminosity density is 
given bv lYounel li 19761) . and we have fitted this density with 
the double power law profile selecting a radial range that 
encloses 99% of the light. The best fit, for 7 and /3 integer 
numbers, gives j — 2, l3 = 5, a = 3.28-Ro, where _Rc is the 
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Figure 1. Fit to the volume density corresponding to a projected 
profile (dashed lines) with the double power law profile of 
eq.|3](solid line). The dotted straight lines in the right panel mark 
the levels 0.5% and 99.5% of the total mass. 



effective radius of the J?^''*-law. Solid lines in Fig. Q show 
this fit, while the dashed grey lines show the i?^/* profile. 

For each gala:xy we have constructed three mass models: 
1) maximum-light (stellar disc + bulge); 2) maximum spher- 
ical DM halo; 3) maximum fiat DM halo. In the first mass 
model, the stellar disc and bulge contribute maximally to 
the rotation curve; in the second and third models, the stel- 
lar mass-to-light ratio is set to the minimum plausible value 
(see Sections l3.ll and l4.lll . The parameters of the dark mat- 
ter haloes are then tuned to reproduce the rotation curve. 
The eccentricity is e = in model 2 and e = 0.95 (axis ratio 
q ~ 0.3) in model 3. 



2.2 Initial conditions 

The orbits of the particles are integrated in the {R, z) plane 
and then uniformly distributed in the azimuthal angle (f). 
Thus each orbit is used to follow an ensemble of clouds 
that are launched from locations lying around a circle in 
the plane. The integration is performed in cartesian coordi- 
nates using the fourth-order Runge-Kutta method. At each 
interval St, the positions and the velocities of the particles 
are projected along the line of sight. We used a variable time 
step with upper value 5i = 3 x 10^ yr. 

Initially, the particles are considered to be gyrating 
around circular orbits in the plane z = in such a way 
that they have isotropic random velocities with dispersion 
o"ran ~ 7kms^^. Then the particle receives a kick velocity in 
a direction that is chosen at random from an axisymmetric 
distribution. The magnitude of v has probability distribu- 
tion 
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Figure 2. Distribution of kick velocities for two values of F = 2 
(top) and F = 10 (bottom) in eq. |21 The two distributions are 
Gaussian with /iv = 75kms~^. 



where 9 is the angle between the direction of the kick and 
the normal to the galactic plane. An natural choice for ^ is 



^g(w, 0) = exp 



2hl cos2r ( 



(6) 



where F and the characteristic kick velocity are param- 
eters to be determined. The highest velocities are confined 
to a cone around the normal to the galactic plane, and the 
width of this cone is set by the parameter F. The functional 
shape of ^ determines the thickness of the gas halo and has 
to be constrained by the observations (Section 13.2^ . In or- 
der to avoid numerous low- velocity integrations, we have in- 



cluded a lower cut-off for kick velocities uthr. 



25kmE 



(5) 



Fig. 121 shows the resulting distributions of kick velocities as 
a function of 6 with — 75kms~^ and two choices of F 
(2, top panel and 10, bottom panel). In section we show 
how to constrain the value of F from the data. It turns out 
that high values (F ^ 5) are required leading to a fairly col- 
limated flux. 

Stochastic acceleration processes sometimes produce a 
tail of particles that extends to large v, which is absent if 
^ is taken to be the Gaussian function ©. Therefore we 
examine the possibility that ^ tends to a power law at large 
V. Specifically, we present results for the choice 
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V ^ Ucut 

V > Vcut, 



(7) 



where u* is a characteristic velocity that we identify with 
the velocity dispersion of the ISM, v, — 7kms^^, and the 
slope Ov and the cut-off velocity Vcut are free parameters. 

We have also experimented with forms of P{v) that are 
not functions of t)/ cos^" 0, as is the case in eqs Q and Q. 
We find, for example, that when P is of the form P(v) = 
(,{v)cos^ 9, the data require even larger values of F than 
those just discussed. 

2.3 Outflow rate and halo mass 

In our model, the flow of material escaping the disc is con- 
sidered constant throughout the life of the galaxy. This is 
a reasonable assumption for a normal spiral like the Milky 
Way given the regularity of the star formation history (e.g. 
iRocha- Pinto et alll2000l) . Therefore, at all times, there will 
be a constant number of particles (mass of gas) in the halo. 
Of these, the number launched between R and R + AR and 
with kick velocities between v and v + Av will be: 



Ani,^io{R,v) = 2-KRhout{R,v)T{R,v)ARAv, 



(8) 



where riout {R, v) is number of particles (outflow rate) gen- 
erated per unit time and area at radius R with kick velocity 
■iT, and t{R,v) is the travel time for these particles. We can 
re-write the outflow rate 



nout{R,v) = nout(-R) X P{v), 



(9) 



where P{v) is the distribution in eq.|^normalized such that 
r°° P(v)dv = l. 

Since the causes of gas outflows are, presumably, strong 
stellar winds and supernova explosions, we assume that the 
outflow rate is proportional to the star formation rate (SFR). 
The latter can be taken proportional to the gas surface den- 
sity Sg(_R) to a power qsf (Schmidt law) and therefore: 

nout(i?) [Eg(i^)]"s^ (10) 

where qsf ~1.3 ^Kennicuttlll989^ and H is a. constant that 
ultimately represents the efficiency of the transfer of kinetic 
energy from the stellar activity to the halo gas. We can re- 
late this parameter to the total mass of the halo gas. If we 
integrate eq.|Hlover the whole galaxy, we find that the total 
mass of the gaseous halo is 



2nmpH J^^''' dR i?Eg J^^ dvT{R, v)f{v) , 



(11) 



where rUp is the mass of the single particle and Rcut is the 
cut-off radius of the SF activity in the disc. Therefore, for 
any given Afhaio we obtain a value for rUpH and an outflow 
rate fea. 110^ in proper units. 

2.4 Comparison to the observations 

The strategy of this study is to compare our models with H I 
observations of spiral galaxies showing the presence of extra- 
planar gas. In the model, at each time interval dt, positions 
and velocities of the particles are projected along the line 
of sight for a specified inclination of the galaxy. From the 



projected velocities we construct a model cube of particles 
that can be compared with the H I data cube. The total flux 
in the model cube is normalized to the total H I mass of the 
galaxy. The spatial and velocity resolutions are smoothed to 
those of the data. A disc component is added to the output 
of the model with a density profile derived from the data 
and a mass Mdiac = Mtot — Mhaio, where Mtot is the total 
H I mass of the galaxy. The number of particles used in the 
simulation is on average ~ 100 per beam. This, in the case of 
NGC891, implies a typical particle mass of ~ 10* Mq, two 
orders of magnitude below the observational detection limit. 
Therefore, our model is smooth and different simulations do 
not show significant fluctuations with respect to each other. 

The surface density Sg(-R) in eg. 1101 has been derived, 
for NGC891, considering both the atomic and molecular 
gas. The distribution of neutral gas was derived from the 
HI data whilst that o f H2 was taken from CO observations 
JSofue fc Nakailll993ll . The flt of the total gas density with 
an exponential profile gave _Rgas=5.0 ± 0.5 kpc and the to- 
tal mass Mhi+H2 ~ 7.9 x 10^ Mq, in good agreement with 
the total mass that one gets using the HI and CO fluxes: 
Mtot = Mm+MH2 = 4.0+(3.2-4.1)x M© , the latter values 
referri ng to two different ch oices for the CO/H2 conversion 
factor isofue fc Nakailll993^ . For NGC 2403, we have con- 
sider only atomic gas, being the dominant component and 
used the exponential fit to the HI surface density given in 
iFraternali ct al. ( 200 



3 APPLICATION TO NGC 891 

The (Sb) spiral galaxy NGC 891 is one of the best known 
and extensively studied nea rby edge-on galaxies. We as - 
sumed a distance of 9.5 Mpc ivan der Kruit fc Searlelll98lf) . 
which leads to a luminosity Lb — 2.6 x 10^" Lq. The disc 
of NGC 891 s hows intensive star f ormation at a rate SFR ~ 
3.8 Mq yr~^ jPopescu et al .1200411 . The halo region has been 
studied at various wavelengths an d shows a variety of com- 
ponents from radio continuum (Allen. Sancisi & Baldwin! 
^7S) to hot X-ray emitting diffuse gas ( Bregman & Pi ldia 
[199^. NGC 891 is also considere d to be very similar to the 
Milky Way Jvan der Kruitlll984^ . 

NGC 891 has been studied in HI several times 
with ever increasing sensitivity (e.g. 'Sancisi & Allen] 
I1979I: lRupenlll99ll: Is^aters. Sancisi fc van der Hulst 19971) 
and most recently with very deep observations (with 
more than 200 hours of integration) obtained with 
the Westerbork Synthesis Radio Telescope (WSRT) 
iOosterloo. Fraternali fc Sancisil l2005l) . We use here these 
new data at 28" (1.3 kpc) resolution. The total HI map, 
rotated 67° counter-clockwise, is shown in Fig. 01 (upper 
left panel). The distribution of neutral gas in the disc of 
NGC 891 is not symmetric, being more extended in the 
South- West side (on the right of the rotated image). The 
location of this extension is not known (it is probably not 
in the line of nodes), so in deriving the rotation curve and 
studying the kinematics, it is often neglected. These new 
WSRT observations show very clearly the presence of ex- 
tended extra-planar emiss ion, much more extended that in 



the previous observations 
;],997). As shown both by 



il997i) and IFraternali et a 



Swaters. Sancisi fc van der Hulsl 



Swaters. Sancisi fc van der Hulsl 



2003), this emission cannot be 
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Figure 3. Two fits to the HI rotation curve of NGC 891 (squares) 
obtained with the maximum fight (top panel) and maximum halo 
(bottom panel) models. The two inner points of the rotation curve 
have not been included in the fit due to uncertainties on their 
interpretation (see text). 

explained by projection of either a flare or a line-of-sight 
warp, but it has to come from a thick component of neutral 
gas surrounding the thin disc and lagging behind it in ro- 
tation. The extent and the lag of the halo gas are the key 
constraints on our model. 

3.1 Mass models for NGC 891 

The rotation curve of NGC 891 (Fig. |^ h as been derived 
using an envelope tracing technique (e.g ISancisi fc AUenl 
11979) applied to the north-west side of the galaxy (details in 
Fraternali, in preparation). The velocity rises very steeply 
in the inner parts of the galaxy {R < 2 kpc), this phe- 
nomenon may be evidence of a fast rotating inner ring or, 
perhaps non-circular motions related to an inner bar (see 
iGarcia-Burillo fc Guelinlll995l) . Due to these uncertainties 
in interpretation of the inner points of the rotation curve we 
have excluded them in our fitting. 

We decomposed the potential of NGC 891 into four 
components: two discs (stellar and gaseous), the poten- 
tials of a bulge and a DM halo. For the stellar compo- 
nent w e have used the photom etry and disc-bulge separa- 
tion of IShaw fc GilmoTq l|l98 9|). w h o use d the same opti- 
cal data as |var^^^^ruitfc Searld (Il98ll) . They found no 



significant difi'erence in accuracy between fits with two or 
three components (the latter with thin + thick discs) . There- 
fore we considered just two components, a relatively thick 
disc (hz — 1 kpc) and a flattened _R^''''-law bulge. This, as 
IShaw fc Gilmorg (^8^ clearly state, is not a perfect de- 
scription of the light distribution in NGC 891 but it is accu- 
rate enough for our purposes. Indeed (see Section 13.411 our 
main conclusions do not depend on the shape of the poten- 
tial. 

Fig. 13 shows the fits to the rotation curve of NGC 891 
obtained with two of the adopted mass models. The first (top 
panel) is a maximum light (disc-|-bulge) model, in which the 
highest possible mass-to-light ratio is assigned to the lumi- 
nosity density. The M/Lb ratios obtained for the stellar disc 
and the bulge are, 3.7 and 6.3, respectively. The disc and 
bulge are dominant all the way out to the last point of the 
rotation curve. The dark-matter halo is pseudo-isothermal, 
having inner and outer slopes of 7 = and P — 2, respec- 
tively (see eq. 

The bottom panel of Fig.j^shows the fit to the rotation 
curve with a maximum-halo model. In this case the M/Lb 
ratios were fixed to 1.5 and 3.5 respectively for disc and 
bulge. We chose these values as the minimum allowed for two 
stellar populations of mean ag es of about 2 and 4-5 Gyr, for 
disc and bulge respectively iCharlot. Worthev fc BressanI 
Il996h . In order to reproduce the rotation curve in this case 
we need to use a more peaked DM halo of a NFW type 
(7 =1, (3 =3). The DM halo dominates everywhere out- 
side R = 2 kpc. The model shown in Fig. |3 is for a round 
halo (e — 0). We also constructed a model for a flat halo 
(e — 0.95, q ~ 0.3), where the fit to the rotation curve is es- 
sentially identical to that for a round halo. The parameters 
of the mass models of NGC 891 are reported in Table Q 

3.2 Model constraints 

The most striking evidence for extra-planar gas in NGC 891 
is given by the total HI map (Fig. 2] upper left panel; see also 
iFraternali et alj i2005ft ) and our first step is to reproduce the 
observed vertical distribution. Hence, we have explored the 
parameter space to obtain the right amount of high-latitude 
gas. The parameters of our model that have a critical impact 
are: 1) the functional form of P{v) (eq.|5Jl with its parame- 
ters (/iv, ctv and Ucut); 2) the mass of the halo gas (eg. IllH : 
3) the cut-off radius of star formation (eg. IllH . 

To constrain these parameters we have used an itera- 
tive method. We have constructed total maps with different 
values of these parameters, subtracted them from the total 
HI map of NGC 891 and minimized the residuals in the halo 
region only (at |z| ~ 1 kpc). For each of the three functional 
forms of P{v), the lower three panels in the left column of 
Fig.BJshow the best fits to total maps that can be obtained 
by adjusting the free parameters. The parameters of these 
models are given in Table |21 In each case the cut-off radius 
for SF is -Rcut = 16 kpc and the potential is that of the 
maximum-light model. Hence the models differ only in the 
distribution of kick velocities. 

The upper right panel of Fig. |1| shows the residuals ob- 
tained just mirroring the data about the minor axis and 
subtracting it from the original total HI map. Clearly, no 
axisymmetric model could produce smaller residuals than 
these. The lower three panels of the right column show the 
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Table 1. Mass models for NGC 891 



Model 




(M/L)disc 






(M/L)bulge 


Re 


^bulge 


PO.DM 
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7 


/3 


EDM 






(M0/L0,b) 


(kpc) 


(kpc) 


(M0/L0,b) 


(kpc) 




(M0kpc-3) 


(kpc) 








maximum 


light 


3.7 


4.4 


1.05 


6.3 


3.3 


0.8 
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4.0 
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round DM halo 
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4.4 
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Figure 4. Comparison between the (rotated) total HI map of NGC 891 (upper left panel) lOosterloo. Fraternali Sc Sancis 3l2005l) and 
the results from three dynamical models. The models all have a maximum light potential but differ on the shape of the distribution 
function of the kick velocities (eq. |5J. The right panels show the residuals of the subtraction of the model total maps from the total HI 
map. The upper right panel shows the results of subtracting a minor-axis mirrored H I map from the original data. In the left panels the 
contour levels are: 0.2, 0.45, 1, 2.5, 5, 10, 25, 5OM0 pc~^. In the right panels the grey background indicates the zero level, darker colours 
imply that the measured intensity is greater than that predicted by the model. Contour levels are the same as in the left panels. 



residuals produced by the model whose total HI map is 
shown on the left. In the plane there are systematic dif- 
ferences between data and models; these residuals reflect 
both inhomogeneities in the plane and uncertainties in the 
shape of the radial HI profile. Inevitably, the model does 
not match the outer part of the disc on the right (originally 
South- West) side of the galaxy. 

Our models tend to produce a flaring shape of the total 



map because particles launched from larger radii can reach 
higher distances from the plane. Below the plane (South- 
East in the original orientation), the data show a hint of a 
flare, but the effect is weaker in the data than in most mod- 
els. This flare can be reduced by increasing the parameter 
asF in eq. 1101 The optimum value was found to be close to 
1.3. However, overall the residuals are insensitive to qsf, so 
for NGC 891 we henceforth simply set asp = 1.3. In Section 



A dynamical 



Table 2. Parameters of the maximum light models for NGC891 



Model 


(kms""*") 




Mhalo 

(109 Mq) 


^max 
(km S""*") 


Gaussian 


75 




2.0 


oo 


Power Law 




3.6 


3.0 


oo 


Power Law + v^ut 




2.4 


2.3 


200 



14.41 we vary asF in an attempt to account for gas at very 
peculiar velocities near the centre of NGC 2403. 

In Fig. 1^ we plot density profiles perpendicular to the 
plane, the vertical ordering of the models being as in Fig. 01 
The left-hand profiles are along the minor axis, while the 
right panels are along the parallel line that cuts the major 
axis 14 kpc to the left of the centre in Fig. 0] Circles show 
the data, and the column-density scale is in Mopc"^ The 
dashed lines show the last contour level in the total HI map. 

The best results are obtained with the Gaussian and 
the Power Law (PL)-I-Vcut models. The PL model without 
a velocity cut-off does not produce much extra-planar gas 
above the detection limit, instead it produces a long tail of 
emission below this limit (partially visible in Fig. This 
is why the total map of the PL model appears thinner than 
the others and there are substantial residuals. Note however 
that a lower value of Ov would have produced too much 
material at very high latitude. This shows that, except to 
produce the high latitude spur visible in the total HI map, 
there is no need for a population of clouds with very high 
kick velocities (above 200kms~^) that may arise from an 
un-truncated PL distribution of velocities. Therefore, in the 
following models, we prefer to use the Gaussian distribution, 
which has only one free parameter. The characteristic kick 
velocity required to reproduce the distribution of the halo 
gas is around 75kms~^ (~ 1/3 of the circular velocity in 
NGC891). This number is weU constrained (ilOkms"^): 
for lower hv the extent of the halo is not reproduced. 

It is interesting to compare this characteristic value 
for the kick velocities with the expansion velocities for 
superbubbles that are pred icted by hydrodynan i ical in od- 
els. In the classic paper bv lMac Low fc McCravl (Il988l) . in 
the case of an exponential atmosphere and reasonable val- 
ues for the other parameters, the expansion velocities at 
heights above 300 pc (when the supershell "leaves" the 
plane) are expected to be of iicxp ^ 50kms~^ in good agree- 
ment with our findings. Moreover, observations of verti- 
cal motions in HI gas give very similar values. A well 
known example is the expanding supershe ll found in M 101 
iKamphuis. Sancisi fc van der Hulstl[T99lh for which a line- 
of-sight expansion velocity of ~ 50kms~^ has been derived. 

We now study the effect of the parameter F in eq. |S| 
This parameter affects the direction of the kick velocities. 
The larger F is, the closer the initial trajectories are to the 
normal to the plane (Fig. |5J . Fig. |S| shows channel maps 
(heliocentric velocities are shown in the upper left corner) 
for the data of NGC 891 and for three models with differ- 
ent values of F. The three models have the same parameters 
as the Gaussian model in Table |5| and Fig. |1] The channel 
maps shown here are those at extreme approaching veloci- 
ties. Varying the opening angle produces different patterns 
especially for these channels. As one can see from this figure. 
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Figure 5. Comparison between the vertical density profiles for 
NGC 891 (circles) and the three models (solid lines) shown in 
Fig. Bl (from top: Gaussian, PL, PL-|-Vcut). The left panels are 
taken at the centre of the galaxy and the right panels at about 
14 kpc towards left (North-East in the original orientation). The 
left parts of the plots correspond to the lower parts in the maps 
in Fig. |1] The column densities are in AIq pc~^ and the dashed 
lines show the level of the last contour in the total H I map in 
Fig.H 
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Figure 6. Comparison of three channel maps of NGC 891 (first 
column) and those of three models with different opening angles 
for the kick velocities (parameter F in eq. |SJ. The first column 
shows the data with heliocentric radial velocities reported in the 
upper left corners (Vsya = 528kms~^). The data are consistent 
with values of F ^ 5. 



if the cone of the kick velocities is too wide (low F) , particles 
end up populating channels that are not populated in the 
data. For F ^ 5 the channel maps are similar to those of the 
data. We have fixed F = 8 for all our models. When F ^ 5, 
P{v) with u ~ /iv falls by a factor two as 9 increases from 
to 15°. Thus, the data require a fairly highly coUimated 
upward flow. 
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Figure 7. Comparison between 5 channel maps for NGC 891 iOosterloo. Fraternali &: Sancisil 12005) and those obtained with three 
dynamical models with different distribution of the kick velocities fTablel^. The first column shows the data, heliocentric radial velocities 
are reported in the upper left corners. The channel map in the bottom row is roughly at the systemic velocity (Vaya = 528kms~^). 
Contour levels (for data and models) are: 0.45 {2cr), 1, 2, 5, 10, 20, 50 mjy/bcam. 



3.3 Dynamics of the extra-planar gas 

Extra-planar gas rotates more slowly than the gas in the 
plane. This may be seen in the channel maps. Fig. |7| 
(leftmost column) shows five representative channel maps 
in the approaching (North-East) side of the galaxy, com- 
pared with the three models of Table |5| The emission of 
NGC 891 becomes thinner the more one approaches high 
rotation velocities (first rows), far from the systemic one 
{vsyB = 528 kms^^). Although this effect can be produced 
also by o ther phenomena, especially a line of sight warp, 
ISwaters7 Sancisi & van der Hulst (1997) showed that the 
shape of the channel maps is accurately reproduced only by 
a lagging halo. The presen ce of a mild warp along the li ne of 
sight cannot be excluded iBecauaert fc Combed [T997ft but 
it will play a minor role. In our models we have assumed 
that the galaxy is perfectly edge-on at each radius, so all 
the emission outside the plane is produced by extra-planar 
gas. 

Although the models reproduce the basic shape of the 
channels, there are systematic differences with the data. In 



the data, extra-planar gas is more concentrated to channels 
close to the systemic velocity (bottom row) and is almost ab- 
sent from the extreme channels (two first rows) . Even though 
the model channel maps at Vhci = 300.2 km s~^ are thinner 
than those at Uhoi ~ 357.9 kms"'^, the models show much 
more extra-planar emission than the data in the extreme 
channels. The Power-Law model (third column) has globally 
less extra-planar emission than the others (see Section f3.21 
but the kinematic pattern is the same. Thus regardless of the 
distribution of kick veloci ties, model extra -planar material 
invariably rotates too fast. lFraternali et al.l 12005) have esti- 
mated the vertical gradient in rotation velocity in NGC 891 
to be about 15kms~^kpc~^ (in the range 1.3 < z < 5.2 
kpc). We have applied the same technique to our model 
cube and found a value of ~ 7kms~^kpc~^ (a factor of 
two smaller than in the data). We have also inspected the 
behaviour of v^{R,z) in our model and found a very simi- 
lar value for the rotation velocity gradient ~ 6kms~^kpc~^ 
between i? ~ 5 kpc and Rcut (the cut-off radius of star 
formation). In the inner regions the gradient is more pro- 
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Table 3. Parameters of the DM-dominated models for NGC891 



Model 








(kms^-*-) 




Round DM halo 


70 


1.85 


Flat DM halo 


80 


2.4 



nounced, beyond Rent the gradient become steeper close to 
the plane (z ^ Ikpc) and then shallower afterwards. The ve- 
locity gradient inferred from the observations relates to the 
region 5 kpc ^ R ^ Rcut and has to be compared with the 
value of 6 — 7kms^^kpc^^ derived from the model. In the 
innermost regions, the data also indicate the possibility of a 
larger gradient (Fratcrnali ct al. 2005). 

It seems that clouds that are shot up loose part of their 
angular momentum, rather than conserving it. A possible 
mechanism for this loss is interactions with pre-existing hot 
halo or infalling material from the IGM (see Section |7|l. 

3.4 Impact of the potential 

We now investigate the effects of different galactic potentials 
for NGC891. The second and third rows of Fig.|Hlshow the 
total maps obtained with our best fit models with potentials 
dominated by the dark matter halo, which may be round or 
flat (see Table [TJ. 

The best fit parameters of the models with maximum- 
DM potentials differ slightly from those obtained with the 
maximum- light potential (see Table I^J. In particular, for 
the first model we have = 70kms~^ and Mhaio = 
1.85 X 10® Mq, both lower than the values required for the 
maximum-light model. This finding reflects the fact that in a 
round potential it is easier for particles to leave the plane of 
the disc. The distributions of extra-planar gas and the resid- 
uals obtained with these two models are comparable. They 
are also comparable with those obtained with the maximum 
disc models. 

What about the kinematics? Fig.|^shows channel maps 
equivalent to those shown in Fig.Q These channel maps dif- 
fer little from those obtained with the maximum-light poten- 
tial. In particular, the upper channels (away from systemic 
velocity) are again thicker than those of the data, so the 
halo-dominated models fare no better in producing the re- 
quired lag in the rotation of the extra-planar gas. Thus the 
dynamics of the extra-planar gas does not depend signifi- 
cantly on the potential of the galaxy. 

3.5 Outflow and inflow rates 

In this section we ask the question whether the model is 
physically feasible in terms of mass loss and energy input. 
In particular we look at the fate of the gas pushed up from 
the galactic plane through the lifetime of the galaxy. Fi- 
nally, we calculate the energy input necessary to maintain 
the fountain. 

The good agreement between the model and observed 
total-Hl maps guarantees that the models have the right 
mass of extraplanar gas. However, the mass of the halo gas 
that we obtain from eq. 1111 is unrealistically large: for all 
models we obtain Mhaio — 2 x lO" Mq , which is half of the 
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Figure 9. Comparison between five obs erved channel maps of 
NGC891 lOosterloo. Fraternali fc Sancisil 120051 and those pro- 
duced with the two dynamical models in Table 1^ The first col- 
umn shows the data, heliocentric radial velocities are reported 
in the upper left corner. The channel maps of the bottom row 
are roughly at the systemic velocity {Vsya = 528kms~^). Con- 
tour levels (for data and models) are; 0.45 {2a), 1, 2, 5, 10, 20, 
50 mjy/beam. 




R (kpc) Vi^i;,^ (km/s) 

Figure 10. Travel times for particles in the maximum light po- 
tential of NGC 891 as a function of initial radius (left) and mod- 
ulus of the kick velocity (right). The dashed line in the left panel 
shows the revolution time in the disc as a function of radius. 



total HI mass of the galaxy and about 4 times larger than 
the "observed" halo mass. This discrepancy arises because 
ea. llll counts everything that is shot up, while the observed 
halo mass includes only material at distances z above the 
plane that exceed the observational resolution, ~ 1.3 kpc in 
the case of NGC 891. Indeed if we consider only the material 
travelling above 1.3 kpc from the plane we obtain a mass of 

Afhalo,z>1.3kpc ~ 6 X 10** Mq. 

It is interesting at each radius to compare the rate at 
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Figure 8. Comparison between the (rotated) total H I map of NGC 891 (upper left panel) iOosterloo. Fraternali fc Sancisill200.'jl and 
those of the two dynamical models in Table O The right panels show the residuals of the subtraction of the model total maps from the 
total H I map. The upper right panel shows the results of subtracting a minor-axis mirrored H I map from the original. In the left panels 
the contour levels are: 0.2, 0.45, 1, 2.5, 5, 10, 25, 5OM0pc~'^. In the right panels the grey background indicates the zero level, darker 
colours imply that the measured intensity is larger than that predicted by the model, and vice-versa for lighter colours. The contour 
spacing is the same as in the left panels. 



Table 4. Inflow/outflow rates and energy input for NGC 891 



Model 


Total outflow 

(A^Qyr-^) 


Net outflow 


Energy input 
(1039 ergs-l) 


Maximum light 








(^>0) 


33.1 


2.1 


49.8 


(z > 1.3) 


5.7 


1.4 


19.2 


Maximum DM 








(^>0) 


22.7 


1.3 


33.2 


(z > 1.3) 


5.6 


1.0 


15.1 



which material is shot up with the rate at which it returns 
to the disc. The amount of mass pushed up is related to 
the outflow rate by eq. |H| Depending on their initial radius 
and kick velocity the particles remain in the halo for a time 
r{R,v). Fig. IIUI shows the travel times for 2000 particles as 
a function of radius (left panel) and kick velocities (right 
panel) for a maximum Ught model. The travel times are 
quite short: half to a third of the dynamical (revolution) 
times at each radius (dashed line). They are slighly larger 
in a rounder (DM dominated) potential. 

Once the mass of the halo has been fixed, the orbits have 
been integrated and the travel times derived, we can infer the 
outflow rates. Fig. 1111 shows the outflow (upper panel) and 
inflow (middle panel) rates at each radius for the maximum- 



light model. Dashed lines show the flow rates at z = 0, while 
solid lines show the flow of material at = 1.3 kpc, and thus 
give an indication of the flow of material that would appear 
as extraplanar in the observations. Especially at small R 
the great majority of particles never reach \z\ = 1.3 kpc, 
and nearly all these particles return to the disc within the 
same radial bin. Near the outside of the star-forming disc, a 
significant fraction of the particles that leave the disc pass 
through l^l — 1.3 kpc. 

The bottom panel of Fig. 1111 shows the rate of flow of 
material across a cylinder of radius i? as a function of R; the 
dashed curve shows the flow for the whole cylinder and the 
full curve shows the flow for \z\ > 1.3 kpc. This flow peaks 
at M ~ 1.4 Mq yr~^ at the edge of the star-forming disc. If 
gas is permitted to cross the plane on its first return to the 
plane, the net fiow is reduced to M ~ 0.2 Mq / yr. We discuss 
such 'second-passage' models in detail in Sections 0] and 13 
Table |1] reports the global values of outflow and inflow rates 
for NGC 891, ai z — and at z = 1.3 kpc. These values are 
very s imilar to those obtained bv lCollins. Benjamin fc Randl 

From the mass outflows we can derive the total energy 
input required to maintain the HI halo, which turns out to 
be a few xl0*°ergs~^. If we consider a supernova rate of 
~ 0.04 yr"^ for NGC 891 (ca lculated from the SFR assum- 
ing a Scalo IMF iScaldll986l) ') and an energy per supernova 
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Figure 11. Rates of outflow (upper panel) and inflow (middle 
panel) per unit area of the disc of NGC891. The bottom panel 
shows the rate at which material flows through a cylinder of radius 
R. Dashed lines include all material, while full lines show only 
material that crosses the surfaces \z\ = 1.3 kpc. The potential is 
for maximum- light. 



of 10^^ erg, the efficiency required for transferring kinetic 
energy to the extra-planar gas would he rj 0.039. There- 
fore the fountain flow required by our model is energetically 
plausible. 

In summary, the application of our model to NGC 891 
has led to the following conclusions: 1) the model repro- 
duces the vertical gas distribution with characteristic kick 
velocities of about 75kms^^; 2) the energy requirement to 
maintain this mechanism is less than 4 percent of the total 
energy input from supernovae; 3) the model's main short- 
coming is that it gives the extra-planar gas too little lag; 4) 
this lack of lag is independent of the choice of potential. 
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Figure 12. Fits to the HI rotation curve of NGC 2403 (squares) 
obtained with our maximum disc (top panel) and maximum halo 
(bottom panel) models. Their parameters are shown in Table ISl 



tori ssen et ai]ll999^ a nd a global SFR rate of 1.2M0yr"^ 
iKennicutt et all200i) . 

NGC 2403 is inclined at 63° rather than being edge-on. 
This intermediate inclination has the disadvantage that we 
cannot directly determine the distance from the plane at 
which emission occurs. It has, however, the advantage that 
it enables us to distinguish radial outflow from inflow. 

Recent H I observations of NGC 2403 have revealed 
a component of neutral gas at lower rotation veloci- 
ties (20 — 50kms~^) with respect to the rotation of the 
thin disc dSchaap. Sancisi fc Swater'3l200d : iFraternali et alJ 
I2OOII) . Three-dimensional modelling of the HI data cube 
has shown that this component is produced by lagging 
extra-planar gas as observed in edge-on galaxies, and re- 
vealed evidence that t he extra-planar gas is flowing inward 
JPraternali et alJl200lll . 



4 APPLICATION TO NGC 2403 

NGC 240 3 is a nearby Sc gala xy located at a distance of 
3.2 Mpc JPreedman et alJl200J) . NGC 2403 is much less lu- 
minous (Lb = 8.2 X 10^ L0) and less massive than NGC 891 
and is in several respects similar to M33 in the Local Group. 
It has a fairly high SFR with very bright HII regions 



4.1 Mass models for NGC 2403 

We have used the H I rotation curve derived by 
IFraternali et alj ||2002|). Since no stellar bulge is detected 
in this galaxv jKen^^SsI) . there are just three contributors 
to the potential: stellar and gaseous discs, and DM halo. 
The scale le ngth of the st ellar and gaseous disc s are 2.0 kpc 
^Kent^ll985^ and 5.7 kpc jFraternaH et alJl2002l) . 
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Table 5. Parameter of the mass models for NGC2403 



Model 


(M/L)di,e 


i?d 




PO.DM 


a 


7 


/3 


EDM 




(Mo/L0,b) 


(kpc) 


(kpc) 


(M0kpc-=^) 


(kpc) 








maximum light 


1.70 


2.0 


0.4 


3.1x10'' 


4.5 





2 





maximum round DM halo 


0.85 


2.0 


0.4 


6.9x10"^ 


14.0 


1 


3 






Table 6. Parameters for the models of NGC2403 



Model 








(kms^-"-) 


(109 M©) 


Maximum disc 


70 


0.5 


Maximum halo 


65 


0.45 


Second passage 


70 


0.5 



Fig. 1121 (top panel) shows the fit to the rotation curve 
obtained with the maximum-disc model for M/Lb ~ 1.7. 
The DM halo in the maximum-disc fit is a double power- 
law model (eq. 3J with 7 = 0, /3 = 2 and a = 4.5 kpc. 
The bottom panel of Fig. 1121 shows the fit to the rotation 
curve of NGC 2403 obtained with our maximum-DM model, 
in which the stellar disc has M/Lb = 0.85, half the value in 
the maximum-disc model. With this choice, the DM halo is 
dominant at all radii. The other parameters of these mass 
models are listed in Table 13 

4.2 Dynamics of the extra-planar gas 

Both displacement of gas from the galactic plane and rota- 
tional lag shift emission towards the systemic velocity and 
away from the velocities of disc gas. This phenomenon is ap- 
parent in a position-velocity (p-v) plot along the major axis 
of the galaxy (Fig. 1131 up per panel). We use here the data 
of iFraternali et alJ l20o3) at a resolution of 30" (0.46 kpc). 
The extra-planar gas in this p-v plot is visible as a tail in 
the region of low rotation velocities towards the systemic 
velocity (indicated by the left arrow). It is difficult to say 
what percentage of this low-rotation gas (said to form a 
"beard" ) is produced b y lag rather than by thickness efi'ects 
jSchaap. Sancisi fc Swa tcrs 2000). There is also halo gas at 
very high velocity, as indicated by the right arrow in Fig. 1131 
This is "forbidden" gas because it shows up in quadrants of 
the p-v diagram (upper right and bottom left) that are in- 
accessible to gas that moves on circles in the same sense as 
the disc, regardless of its distance from the plane or speed 
of rotation. We defer discussion of the nature of this gas to 
Section El 

The lower panels of Fig. 1131 are p-v plots derived with 
three dynamical models. All the models have a Gaussian 
distribution of kick velocities with — 70 or 65kms~^ 
for maximum-disc and maximum-DM, respectively (Table 
|SJl. We have kept these values similar to those found for 
NGC 891 since it is reasonable to think that they are related 
to physical characteristics of the ISM, rather than to the 
potential of the galaxy. 

In Fig. I13l the second panel down shows the maximum- 
disc model, while the third panel down shows the maximum- 
DM model. The mass of halo gas used here is about 
Mhaio — 5 X 10* Mq, considerably lower than that needed for 
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Figure 13. Position-vel ocity plot along the m ajor axis of the 
spiral galaxy NGC 2403 iFraternali et alj|2002l) top panel com- 
pared with the results of three dynamical models. The left ar- 
row in the data p-v plot indicates the location of the extra- 
planar gas, the right arrow points at the forbidden gas (see text). 
In each panel the horizonthal lines shows the systemic velocity 
(Vsys = 133kms-i). Contour levels are; 0.5 (~2.5 <t), 1, 2, 5, 10, 
20, 50 mjy/beam. 

NGC 891. This halo mass is also quite close to that derived 
from the data i.e. about 3 x 10* Mq. Note, however, that the 
technique used to separate the halo gas from the disc is, in 
this case (IFratern ali e^lj|200^. totally different from that 
used for NGC 891 dSwaters. Sancisi fc van der Hulstlll9^7^ 
and a direct comparison may be misleading. 

Fig. 1131 shows that both models reproduce the global 
shape of the p-v diagram along the major axis, although 
they cannot reproduce the forbidden gas in the central parts 
of the galaxy (right arrow in top panel) ; the models produce 
tails only in the region of low rotation velocities, and the 
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Table 7. Inflow/outflow rates and energy input for NGC2403 



Maximum disc 



Second passage 



Model 


Total outflow 

(A^eyr-i) 


Net outflow 
(Moyr-l) 


Energy input 
(10^^ ergs-i) 


Maximum disc 


3.15 


0.23 


3.85 


Maximum halo 


2.72 


0.21 


2.93 


Second passage 


1.70 


0.16 


2.05 



brightness of such tails is very similar to that in the data. 
We stress that in the case of NGC 2403 we did not perform 
any minimization of residuals but simply used values for 
hv close to those obtained for NGC 891 and tuned the halo 
mass and the outer radius of star formation {Rcut — 15 kpc) 
to match the p-v diagram along the major axis. We have 
experimented with different values of these parameters and 
found that is not tightly constrained: one can get a good 
fit to the p-v plot along the major axis with values as low 
as 40kms~^. 

Fig. 1141 compares five representative channel maps for 
NGC 2403 with those obtained with our models. The ar- 
rows in the data indicate features of the extra-planar gas 
that should be compared with the models. The typical lo- 
cation of the extra-planar gas is that shown by the arrow 
in the second row, which is for Vhei = 52.5 kms~^, the ve- 
locity of material that is rotating a shade more slowly than 
the disc gas. The corresponding model maps reproduce this 
emission quite well. The arrows in the bottom two channels 
show peculiar features in the extra-planar gas: in the fourth 
row (uhei = 135.0) a massive 8kpc-long filament, and in the 
fifth row {vhci = 181.4) forbidden gas (see Section f4.4|l . Al- 
though our smooth, axisymmetric models cannot reproduce 
these features, in the third row (whei = 104.1 kms~^) the 
maximum-disc model does reproduce a similar feature. 

The travel times for particles in the potential of 
NGC 2403 are generally longer than those found for 
NGC 891 due to the lower attraction from the disc. As we 
have seen, the mass of halo gas required to reproduce the 
data is smaller. This implies that the outfiow rates are much 
reduced in this galaxy with respect to NGC 891. 

Table |3 shows these values for the three models pre- 
sented here together with the energy input required to main- 
tain the H 1 halo. These values are about an order of magni- 
tude lower than those obtaine d for NGC 891 (Table IH. For 
a supernova rate of 0.01 yr~^ iMatonick et al.lll997l) and an 
energy per supernova of 10^^ erg, the required efficiency is 
r] ^ 0.01. Hence, only a small fraction of the energy released 
by supernovae is required to maintain the fountain. 



4.3 Inflow or outflow? 

We now turn to the most significant discrepancy between 
the data and the models of NGC 2403. In the third row 
of Fig. [H 

an arrow emphasises that at velocities close to 
the systemic velocity (usys = 133kms~^) the emission from 
extra-planar gas is rotated slightly counter-clockwise in the 
channel map. This phenomenon, can be better appreciated 
in p-v plots perpendicular to the major axis. Fig. 1151 shows 
such plots for three perpendiculars: the middle row is for 
the minor axis, while the plots above and below are for 2' 
north-west and south-east of the minor axis, respectively. If 
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Figure 15. Comparison between position-velocity plots per- 
pendicular to the maj or axis of NGC 2403 in the first column 
iiFraternali et all 2002) and two dynamical models. The parame- 
ters of the model are shown in Table The arrows outline the 
asymmetries in the data that may be indication of a general in- 
flow of the extra-planar gas. In each panel horizonthal lines show 
the systemic velocity (Vsys = 133kms~^). Contour levels are: 
0.5 (~2.5 cr), 1, 2, 5, 10, 20, 50 mjy/beam. 



the gas were in circular rotation, these plots would be sym- 
metric with respect to the central vertical line (the location 
of the major axis). Radial motion of the gas inwards or out- 
wards will break this symmetry, as is most easily understood 
by considering the p-v plot along the minor axis. Actually 
all three data plots are strongly asymmetric, as indicated 
by the arrows. To determine whether the observed asym- 
metry corresponds to inflow or outflow, we need to know 
which side of the galaxy is nearer. Absorbing features in the 
disc of NGC 2 403 and the assumption that the spiral arms 
are trailing led iFraternali et alj i200ll) to conclude that the 
south-west side (right in these plots) is nearer. Then the 
asymmetry of the observed p-v diagrams implies that the 
extra-planar gas is flo wing in towards the galactic centre 
jFraternah et alJl200jl . Similar evidence for an inflow of 
the extra-planar has recently been reported for NGC 4559 
( Barbicri ct al. 2005). 

The model p-v plots in the middle column of Fig. lldl are 
slightly asymmetric in the opposite sense to the data plots. 
This phenomenon reflects the fact that in the models there 
is a net outflow of halo gas. It is easy to understand why 
the fountain model predicts an outflow: clouds that form 
the halo are shot upwards nearly perpendicular to the disc 
from platforms that are moving close to the local circular 
speed. Consequently, the clouds start from pericentre. They 
return to the disc after half a period of vertical oscillation. 
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Figure 14. Representative channel maps for NGC 2403 in the first column jFraternali et al1l2002l) compared with the three dynamical 
models in Table 151 Heliocentric radial velocities are given in the upper left corners of each row. The arrows indicate various locations 
of the extra-planar gas (see text). The channel maps in the second row up are roughly at the systemic velocity (Vsys = 133kms~^). 
Contour levels are: 0.5 (~2.5 a), 1, 2, 5, 10, 20, 50 mjy/beam. 



which is comparable to half a period of radial oscillation. 
Hence clouds are near apocentre when they return to the 
disc, and have spent most of the intervening period mov- 
ing outwards. Fig. 1161 illustrates this reasoning by showing 
the {R, z) coordinates along three representative orbits in 
the maximum-disc potential of NGC 2403. Each trajectory 
starts at i? ~ 5 kpc, and from top to bottom the initial ver- 
tical velocity is ~ 70, 140 and 70kms~^. It is striking 
that the outflow inherent in the fountain model produces an 
asymmetry in the p-v diagrams of Fig. ll5l that is appreciably 
smaller than the opposite asymmetry apparent in the data. 
Thus the data imply really significant inflow. 

Is it possible to reconcile the fountain model with the 
observed inflow of halo gas? 

We have seen that ejected gas is bound to flow outwards 
for most of the time prior to its return to the disc plane. 
If the gas disc has a sufficiently low surface density at the 
(large) radius of this impact, a cloud may pass clean through 
the disc. In this case it will contribute to halo-gas emission 
with 7? < for most of the time until it hits the high-surface- 



density inner gas disc. The panels labelled 'second passage' 
in Figs Isl and 14. 21 to 1 1 51 show what the data would look like 
in this case. In particular, the column on the extreme right 
of Fig. ll5l shows that permitting clouds to pass through the 
outer disc restores symmetry to the p-v diagrams for perpen- 
diculars to the major axis of NGC 2403, rather than intro- 
ducing the required sense of asymmetry. This makes perfect 
sense physically, and suggests that our only hope of account- 
ing for the observed inflow with the fountain model involves 
both allowing clouds to pass through the disc, and arguing 
that they are visible only during the second halves of their 
orbits. This condition would be satisfied if outflowing clouds 
were highly ionized (perhaps because they are hot), and the 
clouds become predominantly neutral near apocentre, per- 
haps as a result of radiative cooling (Section I^J. 

4.4 Halo substructures and forbidden gas 

In previous sections we have compared the data of NGC 2403 
with smooth axisymmetric pseudo-data cubes produced 
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Figure 16. Orbits in tlie (R, z) plane for tliree representative 
orbits in tlie maximum-disc potential of NGC 2403. In the bottom 
panel the orbit passes through the disc at the first passage. Right 
panels: solid line=Vz, long dashed line=t)i^, short dashed lme=v^. 



with our models. Although there is a general agreement be- 
tween the two, Fig. 1141 demonstrates that the data contain 
features that are not reproduced by the models. In particu- 
lar several gas complexes are detected as discrete features or 
substructures with masses from about 10^ Mq (limit of the 
data) to 10^ Mq (filament visible in channels at Vhci = 104.1 
and Vhci ~ 135.0 kms~^). Moreover some of these extra- 
planar gas structures are found at anomalous velocities (for- 
bidden gas, see the map Vhei = 181.4 kms~^ and Fig. I13II . 
Can these features be produced in a galactic fountain scheme 
or are they evidence of gas accretion? 

We can think of two ways to obtain these features in 
our models. The first is to increase the kick velocities in 
the centre of the galaxy. This will produce gas with large 
vertical velocities that may end up in the forbidden quad- 
rants of the p-v plot along the major axis. The second is 
to increase substantially the number of particles pushed up 
in the central regions by increasing the parameters asF in 
eq. 1101 Indeed by inspecting the model cubes we found that 
at low density levels some gas is already present at forbid- 
den velocities but below the threshold of the lowest contour. 
Both proposals have reasonable physical explanations. First, 
if many supernovae go off in a limited spatial region (a bright 
HII region), then the supernovae may power a supershell at 
higher expansion velocities. The second idea is also plausi- 
ble considering the uncertainties in the Schmidt law and the 
fact that star formation may well depend on several other 
parameters in addition to gas density. 

The left column of Fig. 1171 shows, from top to bottom, 
two observed channel maps and the major-axis p-v plot. The 
other two columns show the corresponding predictions for 



Dtjla High V, , , High a 




Figure 17. Comparison between two channel maps and a p-v 
plot along the major axis for NGC 2403 (left column) and two 
models. The two models are obtained by increasing the kick ve- 
locity in the centre of the galaxy (second column) and increasing 
the density of fountain particles in the centre (third column) . The 
arrow indicate interesting features in the data (see text). Contour 
levels are: 0.5 (~2.5 a), 1, 2, 5, 10, 20, 50 mjy/beam. 

two models. In the case of the centre column, the ejection 
velocity is = 120kms~^ at i? < 4kpc, comparable to 
the circular speed, Vdvc ~ 130 km s^^. The right-hand col- 
umn shows the result of assuming a high SFR, asF = 3, at 
-R < 4 kpc. In both cases we have enhanced substructure by 
reducing the number of clouds employed in the simulations. 
The top row of panels shows that both models produce fea- 
tures at anomalous velocities that are not unlike some of 
those observed. The second row of panels shows that the 
models can even produce emission at forbidden velocities, al- 
though not to the extent observed (arrow in left-hand panel). 
Unfortunately, the major-axis p-v diagrams of the bottom 
row reveals a serious problem with the models: the data 
show very little gas at speeds larger than circular (arrow at 
left) whereas the models do show such emission, especially 
in the case of large /iv (middle column). 

We summary that there is little promise of accounting 
for very anomalous emission (forbidden gas) within the con- 
text of the fountain model. 

In conclusion, we have shown that 1) the basic distribu- 
tion and kinematics of extra-planar gas in NGC 2403 can be 
reproduced with the same model used for NGC 891 and the 
same characteristic kick velocity (~ 70kms~^); 2) the en- 
ergy input required by this model is less than 1 percent the 
input from supernovae; 3) the models fail to reproduce the 
global infiow of extra-planar gas seen in the data, predicting 
instead an outflow (or no flow); 4) high velocity and mas- 
sive substructures are not easily produced by the fountain 
model. 



5 PHASE-CHANGE MODELS 

The models presented so far have two major failures in re- 
producing the dynamics of the extra-planar gas. They do not 
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Figure 18. Comparison between tliree channel maps for NGC891 in the left column iOosterloo. Fraternali fc SancisJl200?ll and three 
phase-change models. The first column shows the data, heliocentric radial velocities are reported in the upper left corner (Vsys = 
528kms~^). Contour levels (for data and models) are: 0.45 {2a), 1, 2, 5, 10, 20, 50mjy/beam. 



produce enough lag (for NGC 891) and they produce a gen- 
eral radial outflow of the extra-planar gas contrary to what 
is observed in NGC 2403. If these failures cannot be solved 
in a galactic fountain scheme, mechanisms such as accretion 
of material from the IGM, will be required to explain the 
kinematics of the extra-planar gas. Here we explore the only 
remaining possibility to reconcile the fountain models with 
the data. This is that clouds are visible only during the sec- 
ond (infalling) parts of their orbits, because they leave the 
plane highly ionized and only later become predominantly 
neutral. 

We have considered four models of this kind. In the first 
class of models (vertical infall) the particles are switched on 
(become visible as HI clouds) when i = and they at- 
tain maximum height. In this scheme we considered two 
subclasses, with particles stopped at the first or the sec- 
ond passage through the disc. The other class of models 
(radial infall) is constructed by switching on the particles 
at apocentre (R — 0). This conjecture ensures the maxi- 
mum contribution possible from radially infalling material. 
However, only particles with high initial kick velocities reach 
their apocentre before the first passage through the disc (see 
Figure \Wl . Therefore the only radial-infall model that pro- 
duces reasonable results is the one with particles stopped at 
the second passage. 

In Fig. llBl we show the comparison between three chan- 
nel maps for NGC 891 and those obtained with the three 
phase-change models described above. The parameters of 
these models were derived by minimizing the residuals in 
the total HI map and they are reported in Table|S] The cut 
off radius for star formation is i?cut = 13.5 kpc for all these 
models. All the models have maximum-light potentials. 

From Fig. 1181 it is clear that the radial-inflow model 
(rightmost column) is immediately ruled out because it pro- 



Table 8. Parameters of the phase-change models for NGC 891 



Model /iv Mhaio 

(kms-l) (IO^Mq) 



Vertical inflow 75 2.7 

Vertical inflow -|- second passage 80 2.5 

Radial inflow -|- second passage 90 2.4 



Table 9. Parameters of the phase-change models for NGC 2403 



Model hv Mhaio 

(kms-l) (10^ Mo) 



Vertical inflow 70 0.5 

Vertical inflow -f second passage 70 0.5 

Radial inflow -|- second passage 80 0.6 



duces fast rotating extra-planar gas. This is due to the fact 
that beyond the apocentre the azimuthal component of the 
velocity rises again (see bottom panel of Fig. I16|l . The two 
vertical-infiow models are more acceptable. They produce 
slighly more lagging material than the models without a 
phase change (Figs and EJ. This arises because we have 
removed the contribution from the first (outflowing) part of 
the orbits, on which the lag is very limited (Fig. Ilfcil . 

Consider now NGC 2403. Fig. ll9l shows the comparison 
between three channel maps and the same three models used 
for NGC 891. The parameters of these models are shown in 
Table |U] The cut-off radius is J?cut = 11 kpc, and all models 
use the maximum-disc potential. 

The arrows in Fig. 1191 indicate peculiar locations for 
extra-planar gas in the data. In particular the arrow in the 
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Figure 19. Comparison between three channel maps for NGC2403 (left column) and three phase-change models. The first column 
shows the data, with heliocentric radial velocities in the upper left corner (Vaya = 133kms~^). Contour levels are: 0.5 (~2.5 cr), 1, 2, 5, 
10, 20, 50 mjy/beam. 
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Figure 20. Comparison between p-v plots parallel to the mi- 
nor axis of NGC 2403 (left column) and those obtained with two 
phase-change models. The arrows outline the asymmetries in the 
data that may be indication of a general inflow of the extra-planar 
gas. Contour levels are: 0.5 (~2.5 a), 1, 2, 5, 10, 20, 50 mJy/beam. 



bottom channel map (104.1 km s^^) indicates a feature pro- 
duced by radial inflow motions. Clearly only the models with 
particles stopped at the second passage produce a similar 
pattern, whilst the vertical-inflow model displays the oppo- 
site pattern. The same eflFect is indicated by the arrow in 
the channel map in the middle row (52.5 km s"'^) and the 
second-passage models seem to show the correct pattern. 
However, these second-passage models are problematic in 
the top channel map (16.5 km s"'^) in that they show ex- 
tended tails of emission not seen in the data. This emission 
is produce by vertical motions (inflow and outflow) com- 
bined with thickness and lag. This seems to indicate that, 
although the data require a large amount of inflow, the dom- 
inant component of this inflow should be radial rather than 
vertical. 

Finally, consider Fie. 1201 This is the analogue of Fig. 1151 
for the two second-passage models described here. The 
asymmetric pattern shown in the data (evidence of radial 
inflow) is now visible in the model p-v plots. In particular 
in the radial-inflow model which is indeed the most extreme 
inflow model that we can produce. However even in this ex- 
treme (probably unrealistic) model, there is less asymmetry 
than in the data. This shows that the amount of inflow re- 
quired by the data is very large. 

In conclusion, combining the results for NGC 891 and 
NGC 2403, the only promising model is the one with vertical- 
inflow and a second passage because the radial-inflow model 
produces fast rotating extra-planar gas and the model with- 
out second passage shows outflow instead of inflow. However, 
even this vertical-inflow model does not produce as much lag 
(for NGC 891) and inflow (NGC 2403) as the data require. 
Hence, we conclude that selecting portions of the particle 
orbits does not make it possible to reproduce the data with 
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the fountain model. There is an intrinsic need for low an- and this procedure assures a full control of projection and 
gular momentum and radially infalling material that is not resolution effects, 
present in a pure galactic-fountain scheme. 



6 COMPARISON WITH PREVIOUS WORK 

ICorbeUi fc Salpeteillll98i) : ICharlton fc SalpeteJ sug- 
gested that the extended If I discs of galaxies might be built 
up from gas that is shot out from interior, star-forming re- 
gions and later settles to circular orbits of large radius. This 
proposal requires some ejected gas to acquire significant ex- 
tra angular momentum per unit mass about the galaxy's 
symmetry axis, Jz, as it is expelled from the star-forming 
region. We have found that the observations require the out- 
flow to be narrowly coUimated upwards, and in this case 
none of the ejected gas will increase its to the large val- 
ues required for formation of an extended HI disc. Hence 
our work effectively rules out this proposal. 

Our phase-change models are essentially varia nts of the 

BregmanI ljl98nil . In that 
model, star-formation in the disc drives a hot wind that 
cools as it rises, with the result that clouds of H I eventually 
form in it, and fall back to the disc. The main difference 
between our model and Bregman's is that in the latter a 
pressure gradient in the coronal gas steadily accelerates the 
gas over a significant distance, while in our calculations the 
gas is abruptly accelerated within the disc and subsequently 
coasts under gravity. Crucially, in both models the specific 
angular momentum Jz of gas is conserved in the flow because 
the dynamics is axisymmetric. Consequently, in both models 
the rotation speed of gas at cylindrical radius R will be 
Vc{Ro)Ro/R, where Ro is the radius from which the gas was 
shot upwards. It follows that in Bregman's model, as in ours, 
the rotational lag of the halo gas is likely to to be too small. 

Our work close l y resembles that of 
ICollins. Benjamin fc R.andl (|2002h . who also considered 
the dynamics of ballistic gas clouds. A difference in the 
coding is that we integrate orbits in the {R, z) plane and 
then uniformly distribute them in azimuth. This allows 
us to have a larger number of "particles" with a limited 
number of integrations, drastically reducing discreteness 
noise. Howe ver, the main difference b e tween our work 
and that of ICollins. Benjamin fc Randl ll2002h is in the 
comparison with the data. Collins et al. extracted average 
projected velocities from the models, and compared them 
with the mean velocities of ionized gas derived from Ha 
long-slit observations of NGC 5775 and NGC891. The 
problem with this approach is that it does not exploit the 
wealth of information that is contained in the line-of-sight 
velocity distribution. The work of Collins ot al. has been 
recently revised by Jjea ld ct_al, (2 005), wh o applied the 
same ballistic model of ICollins. benjamin fc Randl ll2002l) 
to Fabry- Perot observations of NGC 5775. They found that 
the vertical gradient in rotation velocity predicted by the 
model is shallower than that measured from the data, in 
agreement with our results. 

Finally, we stress that the distinctive feature of the work 
presented here is that we attempt to reproduce directly the 
whole H I data-cube of a galaxy. Our models are constructed 
to match the raw data without performing any data analysis. 



7 CONCLUDING REMARKS 

We have presented a model for the dynamics of the neutral 
gas outside the plane of spiral galaxies (extra-planar gas). 
We have explored the role of internal (stellar activity) pro- 
cesses in the formation and maintenance of the gaseous halo. 
In our model, we have considered orbits of particles expelled 
from the disc by supernova explosions and have integrated 
their paths through the halo. We have generated output 
pseudo-data cubes and compared them with the data of two 
galaxies with a large amount of extra-planar gas (NGC 891 
and NGC 2403). 

The main results are: 

• the vertical gas distribution in NGC 891 is well re- 
produced by the models with characteristic kick velocities 
«kick ^ 70 - 80 km s"^; 

• the same characteristic kick velocity also reproduces ac- 
curately the distribution and lag of the extra-planar gas in 
NGC 2403 (as seen in the p-v plot along the major axis); 

• we were able to constrain the opening half-angle of 
the upward flow for a typical superbubble chimney to be 
e< 15°; 

• the energy input required to maintain the gaseous 
haloes in the two galaxies is < 4 percent of the energy re- 
leased by supernovae; 

• the actual shape of the galactic potential plays a mini- 
mal role in the dynamics of the extra-planar gas; 

• the main problem with our model it that it pro- 
duces only half the required lag of the extra-planar gas in 
NGC 891, so there is a need for low angular momentum ma- 
terial (loss of angular momentum) to reconcile the models 
with the data; 

• a second problem is that our model predicts an out- 
flow of the extra-planar gas whilst, in NGC 2403, the data 
indicate an inflow; 

• high velocity massive substructures, observed in 
NGC 2403 and in NGC 891, are diflicult to obtain in a foun- 
tain model and may be the signature of accretion from IGM. 

• allowing the clouds to be invisible (ionised) for the out- 
flowing part of their orbits does not solve the main problems 
of the models. 

We have seen how our fountain model is able to repro- 
duce accurately some features of the extra-planar gas (such 
as its vertical distribution) and it is energetically consis- 
tent. However it has two main failures that clearly indicate 
that a (collision- less) fountain alone cannot fully account 
for the dynamics of the extra-planar gas. The need for low 
angular momentum material seems to suggest that one has 
to take into account interactions between the fountain flow 
and a hot gaseous halo and/or accretion material. The lat- 
ter possibility seems particularly attractive. The accretion 
of material with low angular momentum (zero on average) 
that interacts with the fountain gas may cause the loss of 
angular momentum necessary to reproduce the lag observed 
in NGC 891. Moreover it may add the component of inflow 
velocity to the extra-planar gas needed to reconcile the pre- 
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dictions of the models with the observed inflow in NGC 2403. 
This possibility will be explored in a forthcoming paper. 

Finally, another topic to pursue is that of our Galaxy. 
The features detected in NGC 891 and other external galax- 
ies are very likely to be the analogous of the Intermediate 
and High Velocity Clouds (IVCs and HVCs) of the Milky 
Way. In the Milky Way we have very detailed information 
about the kinematics of these features, in pa rticular their in- 
flow motion (.Wakker fc van Woerdedll997l) . but also about 
their metallicity that, in t he case of Complex C, has been 
found to be of ~ 0.2 solar dWakker et al.lll999l iTripp et alJ 
|20o3). Such a metallicity, intermediate between a primor- 
dial value and a galactic (roughly solar) value, may indicate 
a mixing process between disc (fountain) gas and the ac- 
cretion material. The study the Milky Way halo will be a 
further application of our model in the near future. 
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